Time reversal symmetry breaking in cuprates 
induced by the spiral spin order. 
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We propose a new interpretation of the spontaneous time reversal symmetry breaking (TRSB) 
observed recently in a pseudogap state of cuprates (Kaminsky et al.). It is shown that the TRSB 
dichroism in ARPES signal may be connected with the local spin spiral structures in system. It 
may be caused by a spin-orbit interaction and by spin polarization of electrons at various sections of 
Fermi surface in spiral state. Angular dependence of dichroism signal is studied in schematic KKR 
approximation. Tests are proposed to check an existence of the local spiral spin structure and to 
distinguish it from the TRSB state with micro-currents constructed by Varma. 
Pacs: 71.10.Fd, 74.20.Rp, 74.20.-z 
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The nature of a pseudogap (PS) state of high-T c 
cuprates in underdoped (UD) region remains the intrigu- 
ing problem {j],^. Recently, using the angular-resolved 
photoemission (ARPES) with circularly polarized light 
(CPL), Kaminsky et al. [|| reveal new property of the 
pseudogap state of UD Bi 2 Sr 2 C'aCu 2 O s s (BSCCO). It 
was shown that this state displays a spontaneous time- 
reversal symmetry breaking (TRSB). Earlier Varma Q] 
had predicted the possibility of TRSB in cuprates. They 
proposed the fascinate ground state with circular mi- 
crocurrents inside the plaquettes of Cu0 2 plane with 
definite alignment of orbital angular momenta associated 
with these micro-currents. Namely, the up-directed or- 
bital momenta arrange along one diagonal and the down- 
directed orbital momenta arrange along other diagonal. 
The proposed in alignment of orbital angular mo- 
menta is not connected with any spin alignment. 

The aim of present paper is to discuss an alternative 
possibility for constructing the state with TRSB. Here we 
propose a state in which the TRSB is due to a spiral spin 
structure. The arguments in favor of such hypothesis are 
following. The electric field of CPL really interacts only 
with the orbital motion. Therefore the TRSB dichroism 
implies a definite orientation of orbital angular momenta 
< L n >^ 0. Such momenta, centered on the atoms, 
may be induced by aligned spin momenta < S n >^ 
through the spin-orbit interaction. This means that in 
ARPES one could observe a TRSB dichroism D if the 
photoemission setup can selectively measure the ejected 
electrons with definite spin projection a =f or J,. Then a 



sign of the TRSB dichroism would depend on a sign of a. 
Since usually the ARPES is non-selective with respect 
to the final spin projection of an ejected electron, the 
summary TRSB dichroism is expected to be zero if a 
mean spin polarization of initial states is zero. But for 
the spiral spin structure the occupancy ri}- a of initial one- 
electron band state {ka} with definite k depends on a. 
Such spin polarization of initial fc-state may induce the 
non-zero TRSB effects in ARPES signal. So we calculate 
a dichroism which might manifest in ARPES in case of 
the spiral spin configuration of cuprates. 

The scheme of experiment [pi is given in Fig. 1. The 
right (left) polarized light with a propagation vector q 7 
impacts the crystal surface determined by a normal vec- 
tor n. The xz plane is one of the mirror planes of crystal 
with x-axis along the CuO bonds or along a diagonal 
direction. The ejected electron has a final momentum 
k/. The ARPES intensity I ~ \M. lf \ 2 S{E t - E f - hu) 
is determined by a matrix element of interaction O = 
(e/2m e c)(Ap + pA) with field 



M = A a F a ; F a =<il) f \p a \il) i {k)> 



(1) 



between initial and final states. In dipole approximation 
it contains the vector potential A of the right (£r = 1) 
or left ((l — — 1) CPL with complex amplitude 

Ar( L ) = A [e x cos9j + iCn(L)ey + e z sin6» 7 ] (2) 

For given configuration of the setup vectors n, k/, q 7 the 
ARPES dichroism signal D is determined by a relative 
difference of intensities for both light polarizations: 
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D = (Mr - M L ) /{M R + M L ) (3) 

We study the symmetry properties of ARPES matrix el- 
ements with respect to reflection in a mirror plane of a 
crystal, which is perpendicular to the surface in a typical 
photoemission experiment. Following ||, consider first 
a time reversal invariant initial state ipi(k) and let q 7 
and n lie in the mirror plane m of crystal (here the xz 
plane). Then, the dichroism signal D is nonzero only if 
kj does not lie in the mirror plane m and D has an op- 
posite signs for k at different sides of mirror plane. Such 
dichroism is called a geometrical one. This large effect 
has been observed at any doping But in UD BSCCO 
the residual dichroism (D ^ 0) has been observed even 
for coplanar configuration of n, k/ , q 7 , in which all three 
vectors lie in mirror plane m. Further we adopt a defi- 
nitions q z ,q and k z , k for the normal and 2D intra-layer 
components of the photon and final electron momenta q 7 
and ky correspondingly. 

Consider first a large geometrical dichroism and then 
discuss a possible origin of the observed residual dichro- 
ism connected with TRSB of the ground state of UD 
cuprate. We suggest that the main contributions to ma- 
trix element give the space regions inside the atomic 
spheres. This is in accordance with a fact that the fre- 
quency dependencies of photoemission intensity roughly 
repeat the dependencies of the photoemission cross- 
sections coming from corresponding atomic components 

The formalism for evaluating the optical matrix ele- 
ment for general lattice within KKR scheme is given in 
0. Some corrections should be introduced to provide 
the common asymptotic ~ e zkfr of the final wave func- 
tion of the ejected electron outside the sample (z > 0). 
We restrict our consideration by the one-step model (see 
||) describing a coherent part of photoemission. Incor- 
poration of rescattering and relaxation processes in frame 
of generalization to a three-step model are needed to de- 
scribe the background in the energy distribution function 
of ejected electron. We believe that one-step model is 
sufficient to describe in qualitative manner the angular 
dependence of dichroism. For this aim we will use the 
most simplified form of the initial and final states in the 



process. 

The starting point in calculating of F a in (1) is the 
KKR wave function for a multicomponent lattice . In 
Hartree-Fock representation the one-particle initial state 
ipi with quasi-momentum k is a superposition of orbitals 
belonging to each center 

<Mr, k) = JL £ B(n z y k ^ l Cl^(r - R n0 ) (4) 

n,/3 

Here /3 enumerates all atoms placed at R n p in unit cell 
n = {fi x , n y ,n z ) and L = l,m are the angular momentum 
quantum numbers of orbitals inside each atomic sphere 
|?' — R n p\ < a<nf3- For the upper valence anti-bonding 
band of CuOi plane the main orbitals ^>\q are the d x 2_ y 2 
orbital of Cu and p x , p y orbitals of two oxygens O x ,O y . 
These orbitals constitute a basis set of the Emery model. 

Thus the initial one-electron state is taken as superpo- 
sition of the main real orbitals and in the second quanti- 
zation representation it is 

^i(k,a) = ^Bi(n,)[c d 4 ff + ic xxL +iCyvL\ ( 5 ) 

n z 

The corresponding site operators d} n(J , x* ncr , y\ a of 
Emery model refer to the real functions d x i _ y i(r — R n ,d), 
p v (r— Rn,v) with Rn,v = -Rn + e^a/2, v = x,y. The func- 
tions are considered to extend inside the corresponding 
atomic (muffin-tin) spheres. The real coefficients c<j, c x , 
c y in (5) are obtained from solution of the Emery model. 
In Eq.(4,5) the term under the sum refer to a layer num- 
ber n z . The amplitudes B(n z ) depending on distance of 
layer from the surface describe in phenomenological man- 
ner a coherent or incoherent interlayer transport along z 
near the surface depending on the phase correlations be- 
tween different layers. For standard bulk initial state tpi 
used in j?J B(n z ) ~ exp (ik z n z ). 

A final state inside the sample is taken in similar KKR 
form with the same in-plane momentum k: 

= E B f (n z y kR ^i l C Lp YZ^ L0 (r - R^) (6) 

Here each of functions ipLp with angular momentum 
quantum numbers L — (I, m) is determined inside atomic 
spheres around corresponding center R n ^. Influence of 
surface at z=0 is described by introducing the factors 
Bf(n z ), by phases 6i,p of complex coefficients 
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C Lp = \C Lp \e s »>, (7) 

and by explicit angular spherical harmonics Yj, = 
Yim(kf) depending on direction of final momentum kf. 
The phases are specific for centers /3 in unit cell and for 
angular momentum /. These phases arise from a match- 
ing of final state (6) inside the sample with the common 
plane wave ~ e lkfT in empty space outside it. The phase 
modulation of contributions in (6) determines the geo- 
metrical dichroism of the photocmission. 

The origin of spherical harmonics in Eq. (6) and of the 
phase modulation of coefficients (7) may be illustrated by 
next consideration. First let us construct the final state 
ipf y p(r) for the electron photoemission into direction kf 
from one center {n, /?} only. According to || it must 
be such function of continuum which has an asymptotic 
form of plane wave ~ e lk/r and incoming radial waves at 
\r — Rn[j\ — ► oo- This final state is 

V>/ = exp (ikfRnp) ^e^Y lm (f)YC m {k f )^(r p ) (8) 

in i 

where rp = \r — R n p\- The scattering phase Sip 
for the orbital momentum I is defined by the asymp- 
totic of real radial function of continuum ipip(r) ~ 
i sin (kr — nl/2 + Sip). (According to the KKR approach 
one may consider that asymptotic form is achieved at the 
surface of the muffin-tin sphere.) 

In similar manner the final state ipf for the electron 
ejected from N s centers of the surface layer should be a 
function which at large z > have an asymptotic form 
with common plane wave ~ e lk/r and incoming spherical 
waves contributed from different centers. 

If kf\R n — R n ',f3'\ > 1 and if we neglect the sec- 
ondary scattering processes, then inside each of nonover- 
lapping muffin-tin spheres surrounding center (n/3) of 
surface layer the final state wave function should just 
have a form (6) with complex coefficients Cl>0 ~ e l8lf) . 
Really the secondary processes synchronize the phases 
5ifj of all contributions into the final state from differ- 
ent angular harmonics and different centers. The KKR 
bulk solution for the final state ipf(kf) found in Q takes 
into account the phase and amplitude synchronization 
of all secondary processes, but it neglect the necessary 
additional synchronization and phase modulation com- 



ing from the boundary surface where solutions should be 
matched with the plane wave with momentum kf . 

In order to study in qualitative manner the angular de- 
pendence of dichroism and its symmetry, it is sufficient 
to use the form (6) of the final state without specifying 
the values and phases of Cl'0 in (6). So we use Eqs.(5,6) 
for schematic representation of initial and final states to 
study the symmetry and possible angular dependence of 
dichroism manifested in ARPES. The components F a of 
the matrix element (1) is expressed as a sum of integrals 
inside atomic spheres around centers (/?) for correspond- 
ing channels / — > /' 

M R(L) = A a (QF^l'k) (9) 

Here A a are components of the vector potential (2) de- 
pending on the right or left polarization of CPL, £ = 
(r(l) — ±1; and k = kf/kf. Functions F^(l', k) corre- 
spond to real initial orbitals v = d x -2_ y 2 1 p Xl p y . To obtain 
them we use the selection rules V — I ± 1 for orbital an- 
gular momenta for integrals inside atomic spheres. For 
simplicity we retain only the matrix elements for the fol- 
lowing transitions p x ( y ) — * s, d and d x 2_ y 2 — > p from 
the O x i y \ and Cu centers of CuOi plane. According the 
KKR calculations such transitions give main contri- 
butions. Omitting of transition d x 2_ y 2 — > f at Cu center 
does not change the symmetry properties of the calcu- 
lated dichroism. This leads only to neglecting the small 
higher harmonics in angular dependence of the ARPES 
intensity. 

The resulting expressions for functions F£(l'k) are pre- 
sented in Table. For the p — ► s, d transitions in oxygen 
they include also the factors 

9x( v ) = s x{y )/^s 2 x + s 2 y ; s x[y) =±sm{k x(y) /2); (10) 

They originate from angular dependence of real ampli- 
tudes c x f y \ and Cd of different orbitals in the initial 
band state (5) of Emery model (with effective param- 
eters ed,e p ,t p d,tpp). At t pp -C t p d the amplitudes in (5) 
are 

c x ( y ) = 9x( v ) shif?; c d = cos?7; (11) 
tan 2ij — 2t p d(cos k x + cos k y )/(ed — e p ). 
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Extension to large t pp does not change the symmetry of 
amplitudes. 

The coefficients C (/j/)(fc) in table include: f) a 
sum over the layers B*f(n z )Bi(n z ) based on the phe- 
nomenological or tight-binding dependencies B(n z ); 2) 
the phase factors exp coming from boundary con- 
ditions; 3) the reduced integrals < V (3 || p \\ 1(3 > over 
angular variables after removing the m dependence; 4) 
the radial integrals; 5) the factors sin 77, cos 77 from am- 
plitudes (11). 

Then the ARPES dichroism signal is 



D{<p) = im{M(AM)*}/(|M| 2 + \AMf 



(12) 



where M = M R + M L , AM = Mr - M L and angles 9, ip 
describe a final momentum . Dependence (12) may 
be presented as D(ip) ~ G(k) sirup with even function 
G(k) relative to reflection in mirror plane zx. According 
to (9) the quantities Mr(j) are determined by complex 
constants Co, Cj, Cjj whereas the rest angular functions 
listed in Table are real functions. It can be shown that 
the dichroism signal is zero if all coefficients Cl/3 in (7) 
have the same phases Spj or their differences are multiple 
to 7r. So a representation of final state with correct phases 
is significant for description of geometrical dichroism. 

At <7 7 = n when photon impacts normally to the C11O2 
plane one obtains 

M = Cu sin 9 cos ip + g x [Co — Cj (cos 2 9 
— sin 2 9 cos 2<p)] + g y Ci sin 2 9 sin 2ip 

AM = -Cjj sin 9 sin ip + g v [C - C/ (cos 2 9 , 

(14) 

+ sin 2 9 cos 2ip)] + g x Ci sin 2 9 sin 2<p 

Here 0,ip are the polar angles of k/. It is seen from 
(11,12) that the dichroism signal D(ip) is odd function 
D{—ip) = —D(p) and it go to zero at ip = n or tp = 0. 
This is an expected property of geometrical dichroism. 

Manifestations of geometrical dichroism depend on nu- 
merous parameters. Fig. 2 illustrates an examples of func- 
tions D(ip) for three angles 9 g — 0, 7r/6, 7t/3 of the 
photon momentum in mirror plane zx and for k mov- 



ing along boundary \k x ± k y \ = 7r and k z 



We 



ascribe the following arbitrary values for relative ampli- 
tudes |Cq/C/| = I Co /Cjj I = 1.0 and the relative phases 



{Sv=o(0), S v=2 (0), S v=1 (Cu)} = {0, 3tt/4, tt/4} of co- 
efficients Co, Cj, Cjj in different channels of oxygen and 
Cu centers. Two setup configurations with x along CuO 
bond or along diagonal direction are considered. Func- 
tion D((p) is odd function of ip and it vanishes at (p = 0, tt. 
The calculated geometrical dichroism disappears for all 
ip if all phase differences Si — 81 = irm are multiple to n. 
This is just the case for the Cu- and O- contributions to 
the matrix element calculated in . There the standard 
KKR bulk wave functions were used and an additional 
phase modulation were neglected. At a normal photon 
impact (0 7 = 0) a dichroism signal is zero on each mirror 
plane of tetragonal lattice, i.e. at ip = \vn. 

Now we take into account the spin-orbit interaction on 
Cu with a constant A 



Vls — A ^ LnS r . 



(15) 



Then the initial band function transforms to ipl+Sip 
in a way equivalent to replacement of d x2 _ y 2 a in (5) by 



id 



(16) 



in Eq.(5,4). Here ^ = a/\a\ = ±1 and C; - X/25E 
where 5E is the energy difference of the d-orbitals of x 2 — 
y 2 and xy, yz, xz symmetries. Additional contribution to 
ip\ a leads to changes M — > M + 5M, AM -> AM + 
8AM(<j) in Eqs.(10, 11,12). The TRSB dichroism signal 
at the normal photon impact ((9 7 = 0) is determined then 
by a quantity 



SAM = £,AC\C U sin 9 cos ip 



(17) 



As a result the dichroism signal D(ip, a) of photoemission 
with the final momentum kf and spin projection a of the 
ejected electron has a form 

D{a, k)=A sin p - ff{^ 2 4C A sin g cos p (18) 

Here M, AM and Cjj are determined by Eqs.(13,14,7) 
and by functions from Table. Only linear in A term is re- 
tained in (16). It is determined only by admixture of d xy 
orbital in Eq.(18). The contributions from d-orbitals of 
xz, yz symmetries in (16) are of a 2-nd order of magni- 
tude in A. The second term in (18) is even function of <p 
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and have nonzero value at ip — or n when all three vec- 
tors q 7 , n, kf lye in the mirror plane xz and geometrical 
dichroism disappears. 

Since a sign of D(ip — 0) depends on sign of spin 
projection a of ejected electron, the overall dichroism 
D = -D(ct, ip) should be zero for initial paramagnet 
(PM) state of system. So, for PM state the dichroism 
at (f = or 7r would be observed only if one selects the 
ejected electrons with definite spin projection on n. For 
this PM state the time reversal symmetry is broken just 
by a measurement of spin polarization of photoelectron. 

However, there are the TRSB states in which the dif- 
ferent regions of k space are characterized by different 
spin polarization. For example, for the ground state 
with a spiral spin structure the TRSB effect manifests 
in ARPES by the nonzero overall dichroism at arrange- 
ment of all vectors q 7 , n, kj in the mirror plane. 

Let us demonstrate the polarization selectivity of the 
level occupancies in fc-space for the spiral state of the 
2D t — t' — U Hubbard model. Calculations were carried 
out for model with U/t = 6, t'/t = 0.1 at doping 0.15 
holes per site. The spiral mean-field (MF) solution is 
characterized by average spins < S„ >= d(e x cos Qn + 
e y sinQn) rotating in the xy plane. We study the MF 
states of two types - with the spirality vectors 

Qi = (it-5Q x ,tt), Qn = (Tr-SQ,TT-SQ) (19) 

along CuO bond or along diagonal direction. The spec- 
tral function A a {k,oj) at to = for definite spin pro- 
jection a on z- axis (perpendicular to the spin rotation 
plane) is equal to 

A <J {k,i) = Y J \<^>f\chMi > ffiEjS^Ei-Ef) (20) 

if 

Here the Fermi function f(Ei) depends on the one- 
electron levels of MF solution and S 7 (x) is the 5-function 
broadened with parameter 7 ~ 0.05i. Fig. 3 presents an 
image of the spin-selective and overall spectral functions 
A a= i(k) and A(k) = J2 a A at cj = for two types of 
spiral states. Dark and light gray lines in Fig. 3 corre- 
spond to main and shadow spin-selective (a =|) sections 
of Fermi surfaces. Similar images for a —I are obtained 
from those for a =\ by inversion k — * — k. The spin de- 
pendence of the level occupancy in k space is connected 



with the spin currents J\ = Jj ~ Q existing in spiral 
state. 

Thus, the TRSB state with spiral spin structure cer- 
tainly has a spin-selective sections of Fermi surface. As 
a consequence, one can observe the TRSB dichroism of 
the ARPES signal even at coplanar arrangement of the 
setup vectors q 7 ,n,kj in mirror plane. Two factors are 
decisive here: 1) The ARPES signal corresponds to a def- 
inite local region of k which for given spirality vector Q is 
associated with definite spin polarization; 2) According 
to Eq.(16,17) a definite spin polarization induces via the 
spin-orbit interaction the orbital angular momenta on Cu 
centers and corresponding nonzero dichroism at coplanar 
setup configuration. 

In order to estimate effect we use the spin-orbit con- 
stant A ~ 800 cm -1 following from excitation spectrum 
of Cu^Cu" and a splitting of d-orbitals SE = 
E x 2_ y 2 — E xy ~ 1 -T- 2 eV in crystal field. Then one has a 
value C\ = 0.0025-^0.005 for amplitude in (16). At setup 
configuration Q 1 = 0, 9 = 8k = 7r/4, (p = tt or and at 
the same arbitrary chosen coefficients (7) as were used in 
Fig.2 flCj/tfol = \Cii/C \ = 1, Sfjj = 0,37r/4,7r/4), one 
obtains for the TRSB dichroism max\D\ = 0.033 -j- 0.066 
at fc corresponding ip = or tt and lying in mirror 
plane. This value is consistent the TRSB dichroism sig- 
nal ~ 3 -j- 5% observed in UD cuprates ||. 

In conclusion, it is shown here that the TRSB dichro- 
ism observed in the ARPES spectra of the UD cuprates 
may be connected with a local spiral spin order in sys- 
tem. This hypothesis differs from the model TRSB state 
proposed by Varma et al. {|,H who connect the TRSB 
with aligned charge circular microcurrents on plaquettes 
of C11O2 plane. Instead, the spiral spin order means 
appearance of the local spin currents J| = — Jj_ of the 
macro-scale, about the domain size. Existence of differ- 
ent domains with different signs and values of the TRSB 
dichroism signal has been shown by study a set of samples 
of UD BSCCO in §. The following test for new hypoth- 
esis may be proposed. The rotation of sample on 180° 
around z-axis does change a sign of the TRSB dichro- 
ism D((p — tt) in our hypothesis and does not change 
the sign in case of the TRSB state constructed in M. 
In the former case the rotation changes the signs of the 
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spin currents and of spin polarization. Relative signs of 
the TRSB dichroism signal at different mirror planes of 
cuprate for two types of spiral states (19) and for the 
state proposed by Varma are illustrated in Fig. 4. These 
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sign B may be measured only if the ARPES signal comes 
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from the same domain of a sample before and after its 
rotation and if the spin currents of spiral state are pined 
during rotation of sample. Note, that in case of the fer- 
romagnet alighnment of spins in the surface layers, the 
TRSB dichroism would have the same signs at all direc- 
tions in mirror planes of the lattice. 

A great sensitivity of the Fermi surface (FS) to the 
spin structure puts the questions important for under- 
standing the pseudogap state of BSCCO: Whether the 
observed FS is a composed result coming from several 
domains with different currents? What is a dynamics of 
these currents and domains in the UD cuprates? May the 
spin fluctuations be frosen near the surface into a static 
domains with spiral or ferromagnet spin order? Addi- 
tional test for the supposed local spiral order is possible. 
One can measure the spin polarization < S > of electrons 
ejected from different sections of the Fermi surface and 
to check the correlations of direction of < S > with the 
sign of the TRSB dichroism signal D(ip) at <p = w. Such 
a program requires a spin-selective detection of photo- 
electrons. 
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TABLE I. Functions F%(l',k) determining the a-com- 
ponents of matrix element in Eq.(l). Index 7 numerates the 
contributions originated from different orbitals of initial state: 



from p x (y\ orbitals of oxygens O 



or cL 



- y 2 1 dxy orbitals 



of Cu; final channels s,p,d correspond to angular momenta 
I' = 0, 1, 2. Functions G x t y \, Lk, g x and g y , are even and 
odd functions relative to mirror planes zx correspondingly. 



Px 
Px 
Py 
Py 
d v 2 



S 
d 

s 

d 
p 
p 



Cog x 
Cig x G x 


CigyPk 

Cu sin cos p 

Cu sin 9 sin p 





Cig x P k 

C Q g y 

Cig y G y 

—Czi sin 8 sin ip 

Cu sin 6 cos p 





Cig x Lk 


Cig y Lh 







*G x (y) = ±sin 2 6>cos2v? - cos 2 + 1/3 
h g x ( y ) are determined by Eqs.(10) 
c Lk — sin 20 cos p\ Pk = sin 2 9 sin 2p 



G 



FIG. 1. The setup configuration of the ARPES experiment 
[3] . The propagation vector q 7 of CPL lies in a mirror plane 
xz and k/, k are the final momentum of ejected electron and 
its component in CuOi plane (xy plane). 

FIG. 2. Dependence of geometrical dichroism D(ip) on the 
azimuthal angle ip of vector k f for k moving along the nesting 
lines \k x ik y \ = n. Solid, dashed and thin curves correspond 
to angles 7 = 0, it/6, tv/3 of the photon momentum q. Setup 
configurations with x along diagonal direction or along the 
CuO bonds refer to left or right graphics. Arbitrary taken 
coefficients (7) are given in text. 

FIG. 3. The images of the spin-selective (plots a,c) and 
overall (plots b,d) spectral functions A-\{k,uj), A(k,u)) at 
u) = in Brillouin zone. Main and shadow Fermi surface are 
shown for the spiral states with the spirality vectors Q — Qi 
or Q — Qn from (19) (plots a,b or c,d correspondingly) . 

FIG. 4. Relative signs of the TRSB dichroism signal for k 
lying on different mirror planes of lattice for spiral states with 
spirality vectors Qi, Qn from (19), or for state proposed by 
Varma (VS). 



